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ABSTRACT 


This  thesis  is  concerned  with  solving  or  approximately  solving  a  maximum-flow 
network-interdiction  problem  denoted  MXFI;  A  network  user  strives  to  maximize  flow 
of  a  commodity  through  a  capacitated  network,  while  an  interdictor,  with  limited  assets, 
attempts  to  destroy  links  in  the  network  to  minimize  that  maximum  flow. 

MXFI  can  be  converted  to  a  binary  integer  program  and  solved  but  this  approach 
can  be  computationally  expensive.  Earlier  work  by  Derbes  (1997)  on  a  Lagrangian- 
relaxation  technique  has  shown  promise  for  solving  the  problem  more  quickly  (Derbes, 
1997).  We  extend  his  technique  and  implement  algorithms  in  C  to  solve  MXFI  for  all 
integer  values  of  total  interdiction  resource  available,  R,  in  some  specified  range; 
interdictable  arcs  require  one  unit  of  resource  to  destroy.  The  basic  procedure  solves 
MXFI  exactly  for  most  values  of  R,  but  “problematic  values”  of  R  do  arise.  For  one  set 
of  test  problems,  a  heuristic  handles  these  values  successfully,  with  optimality  gaps  that 
are  typically  less  than  three  percent. 

We  test  our  algorithms  and  implementations  using  five  test  networks  which  range 
in  size  from  27  nodes  and  86  arcs  to  402  nodes  and  1826  arcs.  Using  a  700  MHz 
Pentium  III  personal  computer,  we  solve  the  largest  problem  in  16  seconds. 
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EXECUTIVE  SUMMARY 


Consider  a  theater  of  war  in  whieh  one  warring  party  uses  a  transportation 
network  to  provide  supplies,  troops,  and  ammunition  to  his  forces  and  suppose  that  the 
opposing  party  has  a  limited  ability  to  attack  and  disrupt  his  opponent’s  use  of  that 
network.  This  thesis  is  concerned  with  solving  or  approximately  solving  a  related 
maximum-flow  network-interdiction  problem  denoted  MXFI:  An  interdictor,  with 
limited  assets,  attempts  to  destroy  links  in  a  network  in  order  to  minimize  the  maximum 
flow  through  that  network  which  an  adversary  might  obtain.  The  problem  is  described  by 
Wood  (1993)  and  was  extended  for  use  at  the  U.S.  Strategic  Command  by  Whiteman 
(1999). 

This  thesis  is  motivated  by  the  possibility  of  weakening  a  military  force  by 
disrupting  its  access  to  supplies,  but  other  uses  may  exist.  For  example,  we  may  wish  to 
disrupt  the  escape  routes  of  a  fugitive  or  reduce  the  flow  of  illegal  drugs  and  precursor 
chemicals  moving  through  a  network  of  rivers,  road  and  air  corridors.  Indeed,  MXFI  was 
originally  motivated  by  the  United  States’  drug  interdiction  efforts  in  South  America. 

MXFI  can  be  converted  to  a  binary  integer  program  and  solved  but  this  approach 
can  be  computationally  expensive.  Furthermore,  the  approach  does  not  lend  itself  to 
efficient  exploration  of  tradeoffs  between  interdiction  resource  expenditure  and 
maximum  post- interdiction  flow.  However,  earlier  work  by  Derbes  (1997)  on  a 
Lagrangian-relaxation  technique  has  shown  promise  for  solving  MXFI  more  quickly  and 
facilitating  the  exploration  of  such  tradeoffs.  We  extend  Derbes’  techniques  in  this 
thesis. 

In  the  Lagrangian-relaxation  technique,  for  fixed  resource  level,  we  relax  the 
interdiction  resource  constraint  in  the  basic  integer  program  using  a  parameter  X.  This 
relaxation  allows  us  to  approximately  solve  the  problem  by  moving  the  resource 
constraint  into  the  objective  function.  The  resulting  problem  is  almost  as  easy  to  solve  as 
if  resource  constraints  were  ignored;  It  is  merely  a  maximum-flow  problem.  We  extend 
this  technique  and  implement  algorithms  in  C  to  solve  MXFI  for  all  integer  values  of 


XV 


total  interdiction  resource  available,  R,  in  some  specified  range.  Interdictable  arcs  are 
assumed  to  require  one  unit  of  resource  to  destroy.  The  basic  procedure  solves  MXFI 
exactly  for  most  values  of  R,  but  “problematic  values”  of  R  do  arise.  A  heuristic  usually 
handles  these  values  successfully:  For  one  set  of  test  problems,  the  heuristic  yields 
relative  optimality  gaps  that  are  typically  less  than  three  percent.  Relative  gaps  are 
typically  small,  a  few  percent,  but  can  be  large  if  the  post-interdiction  maximum  flow  is 
small. 

We  test  our  algorithm,  coded  in  C,  using  five  test  networks  ranging  in  size  from 
27  nodes  and  86  arcs  to  402  nodes  and  1826  arcs.  All  tests  are  performed  on  a  700  MHz 
Pentium  III  personal  computer  with  The  Microsoft  Windows  Millennium  operating 
system  and  Microsoft  Visual  C++  compiler.  The  largest  network  is  solved  in  16  seconds 
for  58  potential  values  of  R. 
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I.  INTRODUCTION 


In  the  maximum-flow  network-interdiction  problem  (MXFI),  an  “interdictor” 
attempts  to  destroy  parts  of  a  capacitated  network  using  limited  interdiction  resources,  so 
as  to  minimize  the  maximum  flow  that  a  “network  user”  can  move  through  the  network. 
This  thesis  investigates  an  efficient  Lagrangian-relaxation  technique  that  solves  MXFI,  at 
least  approximately,  for  all  values  of  a  single  interdiction  resource  in  a  specified  range. 


A,  OVERVIEW 

This  thesis  is  concerned  with  solving  or  approximately  solving  a  maximum-flow 
network-interdiction  problem  described  by  Steinrauf  (1991)  and  Wood  (1993),  and 
extended  for  use  at  the  U.S.  Strategic  Command  by  Whiteman  (1999).  The  problem  is 
defined  on  a  directed  network  G  =  {N,A)  which  has  arc  capacities  Uk  for  all  arcs  k  e  A,  a 
source  node  s,  a  sink  node  t,  and  an  artificial  “return  arc”  a  =  {t,s)  included  in  A  such  that 
Ua  =  °°.  We  are  interested  in  solving  he  problem  for  various  amounts  of  interdiction 
resource  R,  so  we  parameterize  the  problem  by  R: 

MXFI(i?)  z{R)  =  min  max  y 

xeX  y 

s-t.  =0  \fieN 

kGFS{i)  kGRS{i) 

^w^(l-x,)  ykeA 


where  X  = 


{0,1}'''' 


keA 


,  ric  is  the  amount  of  resource  required  to 


interdict  arc  k,  R  is  the  total  amount  of  resource  available,  FS{i)  (forward  star  of  node  i)  is 
the  set  of  arcs  directed  out  of  node  i,  and  RS{i)  (reverse  star  of  node  i)  is  the  set  of  arcs 
directed  into  node  i.  For  fixed  x,  the  inner  maximization  problem  is  a  standard 
maximum- flow  problem  with  arc  capacities  m^(I-x^)  .  But,  the  interdictor  controls  the  x 


variables  and  attempts  to  minimize  that  maximum  flow:  When  Xk=  1,  arc  A:  is  interdicted 
and  its  capacity  goes  to  0;  otherwise,  the  arc  is  left  untouched  and  takes  on  its  normal 


1 


capacity  Uk.  We  consider  only  the  ease  in  whieh  all  are  integer.  Thus,  R  may  be 
assumed  to  be  integer. 

MXFI(i?)  ean  be  eonverted  to  a  binary  integer  program  and  solved  using  standard 
teehniques,  but  this  approaeh  ean  be  eomputationally  expensive  (Wood,  1993).  Earlier 
work  by  Derbes  (1997)  on  a  Lagrangian-relaxation  teehnique  has  shown  promise  for 
effieient  solutions,  but  he  only  solves  MXFI(7?)  for  fixed  R  and  we  believe  planners  will 
often  want  to  explore  the  effects  of  varying  R.  That  is,  they  will  want  to  understand  how 
inereasing  interdietion  effort  results  in  improved  interdietion  results  over  a  range  of 
interdietion  resouree  values.  Therefore,  we  adapt  Derbes’  method  to  solve  MXFI(7?)  over 
a  wide  range  of  values  R.  In  partieular,  we  solve  over  the  widest  range  of  “sensible” 
values,  whieh  are  from  7?  =  0  to  the  smallest  value  off?  that  results  in  a  maximum  flow  of 
zero  through  the  network.  Our  method  typieally  solves  MXFI(7?)  exaetly  for  most  values 
of  R,  but  not  for  a  few  “problematie”  ones.  We  provide  an  effeetive  heuristie  to  deal  with 
these  problematie  instanees. 

B,  LITERATURE  SEARCH 

The  study  of  network  interdiction  began  during  the  Vietnam  War  with  efforts  to 
destroy  enemy  supply  lines  (Wollmer  1964,  1970).  Eater,  the  desire  to  interdict  illicit 
drugs  generated  additional  interest;  see  Steinrauf  (1991)  and  (Phillips  1992). 

MXFI  ean  be  viewed  as  having  evolved  over  time  through  the  studies  of  Wollmer 
(1964,  1970),  Durbin  (1966),  MeMasters  and  Mustin  (1970),  Helmbold  (1971),  Ghare, 
Montgomery,  and  Turner  (1971),  Fubore,  Ratliff  and  Sieilia  (1971,  1975),  Steinrauf 
(1991),  Phillips  (1992)  and  Wood  (1993).  Most  of  these  papers  are  based  on  the 
pioneering  work  on  maximum  flows  by  Ford  and  Fulkerson  (1956). 

Steinrauf  (1991)  solves  the  network- interdietion  problem  with  mathematieal- 

programming  teehniques.  Two  mathematieal  programs  are  developed  whieh  determine 

strategies  to  interdiet  a  network  using  limited  resourees.  The  first  model  identifies  a  set 

of  ares  whose  interdietion  minimizes  the  maximum  flow  through  the  network  while  not 

exeeeding  available  interdietion  resourees.  This  is  essentially  MXFI.  The  seeond  model 

identifies  a  set  of  ares  whose  interdietion  isolates  a  large  set  of  nodes  around  a  speeified 

node,  whieh  might  represent  the  most  likely  loeation  of  a  drug  laboratory.  The  models 
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are  applied  to  a  sample  network  that  is  similar  to  a  river  and  road  network  in  Bolivia 
where  eounter-nareotie  interdiction  operations  were  being  conducted  at  the  time. 

Philips  (1992)  presents  several  pseudo-polynomial  time  algorithms  using 
dynamic  programming  for  interdiction  of  undirected  planar  networks. 

Wood  (1993)  develops  integer-programming  models  for  MXFI  and  variants  and 
shows  MXFI  to  be  NP-complete. 

Reed  (1994)  derives  an  integer-programming  model  to  maximize  the  longest  path 
in  a  PERT  network  through  interdiction.  The  purpose  is  to  slow  the  advance  of  nuclear- 
weapons  technology  in  adversarial  countries. 

Wollmer  (1970)  and  Washburn  and  Wood  (1994)  study  game -theoretic  network 
interdiction  models,  which  are  inappropriate  for  the  problem  we  consider. 

Cormican  (1995)  solves  MXFI  using  Benders  decomposition  and  improves 
computational  speed  with  a  “flow-dispersion  heuristic.”  Cormican  begins  the 
investigation  of  probabilistic  versions  of  MXFI  where  interdiction  successes  or  arc 
capacities  are  uncertain.  Cormican,  Morton  and  Wood  (1996)  develop  this  topic  further 
and  solve  such  problems  with  a  sequential-approximation  algorithm. 

Derbes  (1997)  follows  a  mathematical-programming  approach  for  solving  MXFI. 
He  shows  that  a  technique  based  on  Lagrangian  relaxation  can  be  effective  in 
approximately  solving  the  problem.  It  turns  out  that  the  network  interdictor’s  problem  of 
minimizing  the  maximum  flow  through  the  network  is  difficult  to  solve  because  of  the 
interdiction  budget  constraint.  Therefore,  he  relaxes  this  constraint  using  Lagrangian 
relaxation,  which  allows  the  interdictor  to  violate  the  constraint  while  paying  a  penalty. 
For  a  fixed  value  of  a  penalty  parameter,  the  relaxation  is  an  easy-to-solve  maximum- 
flow  problem  with  a  solution  that  provides  a  lower  bound  on  the  optimal  solution  to  the 
original  problem.  A  corresponding  interdiction  solution  is  also  derived  although  this  may 
or  may  not  be  feasible.  He  maximizes  the  lower  bound  using  binary  search  on  the  value 
of  the  penalty  parameter,  solving  a  maximum  flow  problem  at  each  step  and  guaranteeing 
that  at  least  one  feasible  solution  will  be  found.  The  best  feasible  solution  obtained  in 
this  process  is  taken  as  an  approximate,  possibly  optimal,  solution  to  the  problem,  and  the 
corresponding  upper  bound  is  compared  to  the  maximized  lower  bound  to  judge  solution 
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quality.  Derbes  implements  the  procedure  in  C  and  demonstrates  its  effectiveness.  (He 
also  studies,  with  less  success,  a  dynamic  version  of  MXFI  where  flow  requires  time  to 
pass  through  a  network’s  arcs.) 

Golden  (1978),  Israeli  (1999),  Israeli  and  Wood  (1999)  and  Wevley  (1999)  study 
another  category  of  network  interdiction  problem.  In  this  problem,  an  interdictor 
attempts  to  interdict  (destroy  or  lengthen)  arcs,  using  limited  interdiction  assets,  to 
maximize  the  length  of  a  shortest  s-t  path.  While  these  models  are  similar  in  spirit  to 
MXFI,  that  research  is  not  relevant  to  this  thesis. 

Whiteman  (1999)  presents  a  comprehensive  strategy  to  plan  the  interdiction  of 
complex  infrastructure  networks  that  can  be  modeled  as  capacitated  flow  networks.  The 
strategy  is  designed  to  achieve  a  high  level  of  interdiction  in  a  single  set  of  strikes, 
reducing  the  necessity  of  costly  follow-up  strikes.  He  adapts  the  integer  program  of 
Wood  (1993)  to  select  target  sets  of  minimum  size  and  uses  ad  hoc  Monte  Carlo 
techniques  to  handle  uncertain  interdiction  effects. 

Akgun  (2000)  studies  the  “X-group  network  interdiction  problem”  in  which  a 
network  user  attempts  to  maximize  flow  among  three  or  more  groups  of  nodes,  while  an 
interdictor  interdicts  network  arcs,  using  limited  interdiction  resources,  to  minimize  that 
flow.  This  model  is  more  general  than  MXFI  but  an  extension  of  our  techniques  might 
be  of  value  in  solving  that  problem. 


C.  OUTLINE  OF  THESIS 

The  first  chapter  of  the  thesis  has  introduced  the  max-flow  network-interdiction 
model  (MXFI),  and  given  an  overview  of  earlier  research  in  this  area.  Chapter  II 
provides  definitions,  notation  and  detailed  background  on  MXFI,  and  describes  some 
other  basic  models.  In  Chapter  III,  we  develop  our  Lagrangian-relaxation  approach, 
which  attempts  to  solve  MXFI  for  every  value  of  “reasonable”  R,  and  describe  a  heuristic 
that  is  added  to  deal  with  problematic  values  of  R.  Chapter  IV  provides  computational 
results.  Conclusions  and  recommendations  for  future  work  are  covered  in  Chapter  V. 
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II.  PRELIMINARIES 


This  chapter  provides  essential  definitions  and  notation  for  MXFI(i?),  and 
outlines  the  derivation  of  Derbes’  Lagrangian-relaxation  teehnique  for  solving  MXFI 
(Derbes  1997).  The  notation,  conventions  and  models  follow  Wood  (1993),  Cormican 
(1995)  and  Derbes  (1997). 

A,  DEFINITIONS  AND  NOTATION 

MXFI(7?)  is  defined  on  a  direeted  network  G  =  (N^A)  where  A  is  a  set  of  n  nodes 
and  A  is  a  set  of  m  direeted  arcs.  An  arc  k  =  (ij)  originates  from  tail  node  i  and  terminates 
at  head  node  j.  The  forward  star  of  node  i,  denoted  FS{i),  is  the  set  of  arcs  directed  out  of 
node  i  and  the  reverse  star  of  node  i,  denoted  RS{i),  is  the  set  of  arcs  directed  into  node  i. 

The  network  G  =  {N,A)  has  arc  capacities  W/t  >  0,  a  source  node  5,  a  sink  node  t, 
and  an  artifieial  return  are  a  =  {t,s)  ineluded  in  A  such  that  Ua  =  °°.  An  interdictor  has  a 
total  of  R  units  of  resource  available  for  interdiction. 

A  single  type  of  interdiction  is  assumed  in  this  thesis  and  only  ares  are  interdicted. 
Any  “interdictable  arc”  k  requires  >  0  units  of  resource  to  interdiet;  is  assumed  to  be 
integral  and  thus  R  may  also  be  assumed  to  be  integral.  “Uninterdictable  ares,”  which 
always  include  k  =  a,  have  rk  =  These  ares  may  not  be  interdicted  at  any  eost  for 
political,  tactical,  theoretical  or  other  reasons.  Standard  network  transformations  enable 
the  modeling  of  undireeted  networks  and  node  interdiction  if  desired  (e.g.,  Ahuja  et  al. 
1997,  pp.  38,  46).  Other  generalizations  sueh  as  “partial  interdiction”  (Wood  1993)  are 
possible,  too. 

A  out  {Ns,  Nt),  also  denoted  C,  is  a  partition  of  the  node  set  N  into  two  subsets.  Ns 
and  Nt,  sueh  that  s  e  Ns  and  t  G  Nt.  Eaoh  cut  defines  a  set  of  ares  that  have  one  endpoint 
in  Ns  and  other  endpoint  in  Nt.  With  respect  to  the  out,  an  are  k={i,j)  is  a  forward  are  if  i 
G  Ns  and y  G  W ;  otherwise  it  is  baokward  are.  The  set  of  forward  arcs  is  denoted  by  Ac, 
and  we  often  refer  to  this  set  as  “the  cut’  ’  for  the  sake  of  simplioity.  The  oapacity  of  the 
out  is  .  A  minimum  cut  is  a  cut  whose  oapaoity  is  minimum  among  all  possible 

keAc 
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cuts  in  the  network.  The  maximum- flow  minimum-out  theorem  (Ford  and  Fulkerson 
1956)  states  that  the  maximum  flow  s-t  flow  in  a  network  G  equals  the  capacity  of  a 
minimum  cut. 

B.  MAXIMUM-FLOW  NETWORK  INTERDICTION  MODELS 

1,  Maximum-Flow  Network-Interdiction  Model  Reformulation 

Cormican  (1995)  reformulates  MXFl(i?)  for  easy  conversion  to  a  MIP: 

MXFI-RF(i?): 

Indices: 

ij  e  N  nodes  of  direoted  network  G  =  {N,A)  inoluding  two  speeial  nodes, 
the  souree  s  and  the  sink  t 

ke  A  direeted  ares  in  the  network  G  =  {N,A),  k  =  {ij) 

Data: 

Uk  nominal  eapaeity  of  are  k 

rt  amount  of  resouree  required  to  interdict  arc  k 

R  total  amount  of  resource  available  to  the  network  interdietor 

Decision  Variables: 

Network  user: 

yk  amount  of  flow  on  are  k 

Network  interdietor: 

Xk  1  if  are  k  is  interdicted;  0  otherwise 

Formulation: 

z(i?)  =  minmaxy^-^x,y, 


s.t. 
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V/  G  N 


T^yk  -  Y.yk  =0 

k&FSd)  k&RSd) 


0  <  \/k&  A 


where  X  =  <  x  g  {0, 1}''^'  ^  =  0 


keA 


2.  Integer  Program 

Taking  the  dual  of  the  inner  maximization  yields  the  following  mixed-integer 
program  first  derived  through  other  means  by  Wood  (1993). 


MXFI-IP(i?); 

Indices  and  Data:  As  in  MXFI-RF(i?) 
Decision  Variables: 


Xk  1  if  are  k  is  interdicted;  0  otherwise 

ai  at  =1  if  /  G  Nt,  else  «;■  =  0 


Pk  =  1  if  A:  is  a  forward  arc  of  cut  and  not  interdicted,  else  Pk  =  ^ 

Formulation: 

z(i?)  =  min^M,^, 

"’P’*  ^ 

s.t. 

a.  -a.-Vx^-^  Vk  =  (ij)  e  A -  a 


ct/G  {0,1}  Vi  G  A 

a,.  =  0,  af  =  1 

Pke{0,l}  VkeA 

Pa^O 

Xk  &  {0,1}  \fke  A 


Xa=0 


1 


(1) 


k&A 


The  variables  «;■  identify  a  eut  {Ns,Ni)  with  i  G  Ns  if  cCi  =  0  and  i  g  Nt  if  cCi=  1 .  The 
variables  Xk  and  Pk  represent  interdietion  deeisions  with  respeet  to  {Ns,Nt)  and  have  the 
following  interpretation:  For  forward  ares  k  =  (ij)  in  the  eut,  «’,  -«y=-l  so 

^k+  Pk  is  required.  So,  either  =1,  indieating  that  this  are  is  interdieted,  or 
=1,  indieating  that  this  are  is  not  interdieted  and  forms  part  of  the  minimum  eut  after 
interdietion.  =  0  indieates  that  are  k  is  neither  interdieted  nor  part  of  the 

minimum  cut  after  interdiction. 

3,  Lagrangian  Relaxation  for  the  Integer  Problem 

The  basic  integer  program  formulation  MXFI-IP(i?)  is  hard  to  solve  but  becomes 
easy  if  we  ignore  the  interdiction  budget  constraint  (1).  While  this  constraint  cannot  be 
ignored,  it  can  be  relaxed.  Following  Derbes  (1997),  we  use  Lagrangian  relaxation  and 
moved  the  interdiction  resource  constraint  into  the  objective  using  a  parameter  A.  The 
resulting  problem  is  almost  as  easy  to  solve  as  if  constraint  (1)  were  ignored.  The 
Lagrangian  relaxation  of  MXFI-IP(7?)  is: 


LR(X,R): 

z(/l,7?)  =  miny  +Ar^x^)-JiR 
s.t.  ct,  -aj  +  Xi^  +  Pi^  >0 

a, 

cci&  {0,1} 

as  =  at  =  \ 

Ag  {0,1} 
A=0 

Xk&  {0,1} 


VA:  =  {i,j)  e  A  -  a 


Vi  G  N 


Vke  A 


Vke  A 


Xa  =0 
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It  can  be  shown  that  the  LP  relaxation  of  this  model  has  integer  extreme  points  so 
we  may  take  its  dual  and  simplify  to  obtain: 

LROg/?): 

z{A,  R)  =  max  -  AR 

y 

s-t.  Xk-t  -  Z-fi  =0  V/G  iV 

keFS(i)  keRS(i) 

0  <  <  min{Mj.,ylr^}  VkeA 

LRD(yl,i?)  is  essentially  a  maximum  flow  model  defined  on  G  but  using 
eapaeities  modified  by  the  dual  eost  of  interdieting  the  are. 

Proposition  1:  Any  solution  to  LRD(yl,i?)  finds  a  minimum-out  (Ns,  Nt)  that 
oorresponds  to  a  feasible  or  infeasible  solution  (x,a,P)  to  the  original  problem  MXFI- 
IP(i?)  as  follows: 

1.  a.  =\  Vi  G  Nt,  Oi  =0  V/  G  Ns 

2.  A  =  1  if  i  G  Ns,j  G  N  and  y^.  =u^,  i.e.,  if  A:  is  a  forward  are  of  the  minimum 

out  and  ,  then  k  is  not  interdioted. 

3.  =  1  if  i  G  Ns,j  G  Nt  and  y^  =  Ar,^,  i.e.,  if  A:  is  a  forward  arc  of  the  minimum 
cut  and  u^>  Ar^.  ,  then  k  is  interdioted. 

4.  Xj^  =  =  Q  '^k  that  are  not  forward  ares  in  the  cut.  | 

Note  that  A  oan  always  be  perturbed  so  that  =  Ar^.  does  not  ooour.  The 
solution  X  is  feasible  if  the  interdiotion  budget  oonstraint  (1)  holds. 
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III.  NETWORK  INTERDICTION  BY  LAGRANGIAN 

RELAXATION 


This  chapter  presents  our  approaeh  to  solving  MXFI(i?)  through  LRD(yl,i?),  for  a 
range  of  resouree  values  R. 


A,  OVERVIEW 

We  use  LRD(/l,i?)  as  our  main  model  to  develop  the  algorithms  to  solve 
MXFI(i?).  The  funetion  z{A,R),  whieh  LRD(yl,i?)  eomputes,  is  a  pieeewise-linear  coneave 
funetion  in  A  and  Derbes  (1997)  shows,  for  a  fixed  value  of  R,  how  to  maximize  this 
funetion. 


Beeause  of  the  relaxation,  z{A,R)  <  z{R),  but  for  some  A  values  we  may  find  a 
Lagrangian  multiplier  sueh  that  =  R  ,  where  the  xt  are  eomputed  via  Proposition 


he  A 


1 .  In  this  ease, 

z{A,R)  =  min^M,;^,  +A{Y,V,-R)  =  min^w,^,  =  z{K). 


(2) 


keA  keA  keA 

That  is,  we  have  solved  MXFI(R)  exaetly.  Derbes  (1997)  uses  binary  seareh  on  A  values 
to  maximize  z{A,R)  for  fixed  R  and  often  finds  A  sueh  that  z{A,R)=  z{R),  but  not 
always. 


In  this  ehapter,  we  adapt  Derbes’  method  to  attempt  to  solve  MXFI(7?)  for  all 
integer  values  of  R  over  the  largest  meaningful  range,  whieh  is  from  the  smallest  value 
that  allows  no  positive  post-interdietion  flow  down  to  7?  =  0.  (In  faet,  we  seareh  this 
range  from  large  to  small.)  For  the  sake  of  simplieity,  we  assume  from  hereon  that  rk=  I 
for  all  k;  even  this  speeial  ease  of  MXFI  is  NP-eomplete  (Wood  1993). 


B,  SOLVING  THE  RELAXED  MODEL 

Proposition  2:  Given  a  fixed  A,  with  A  for  any  k,  we  ean  solve  MXFI(7?), 
for  some  as  yet  unspeeified  R,  through  LRD(yl,7?),  as  follows: 
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1.  Find  a  minimum-capacity  cut  Ac  in  G  with  arc  capacities  min{M^,ylr^}  VA:  G  A 

where  =  1  for  all  interdietable  ares  k,  and  rk  =  otherwise.  This  is  easily 
aeeomplished  with  an  effleient  maximum-flow  algorithm,  whieh  eomputes  J[A)  (e.g., 
Ahuja  et  al.  1993,  Chap. 7). 

2.  Ares  in  the  minimum-eapaeity  eut  Ac  are  examined  to  find  the  eorresponding 

“interdietion  set”,  Ac,  where  ke  Aj  if  and  only  if  A:  is  a  forward  are  of  the  minimum 
eut  sueh  that  • 

3.  Let  R  =  \  Aj\  and  define  Xk{R)  =  1  for  all  ke  Aj  and  Xk{R)  =  0,  otherwise.  x{R) 
solves  MXFI(7?);  see  Proposition  1  and  Equation  (2)  above.  We  have,  in  effeet, 
optimized  the  Lagrangian  funetion  z{A,r)  through  LRD(/l,r)  for  R  =  r.% 


So,  we  now  know  that  for  any  fixed  value  of  A  we  ean  obtain  a  solution  to  MXFI 
for  some  value  of  R.  But  we  are  interested  in  solving  MXFI  for  all  (reasonable)  integer 
values  of  R.  Given  how  A.  affeets  the  solution  x(7?)  through  Proposition  1  and  Proposition 
2,  it  is  elear  that  the  only  values  of  A  we  need  be  eoneerned  with  are  in  the  range  {Q,UmaA\- 
Furthermore,  from  Proposition  2,  we  ean  see  that  the  solution  x(7?)  ean  only  ehange  when 
the  definition  of  ^4/  ehanges;  it  ean  be  shown  that  these  are  the  eorner  points  at  whieh  the 
slopes  of  z{A,R)  andy(zl)  change.  Intuitively,  this  is  likely  to  oeeur  at  A=  ut  for  some  k, 
although  examples  ean  be  given  to  show  that  not  all  eorner  points  oeeur  at  sueh  A. 
Nonetheless,  we  ean  solve  many  problems  with  aeeeptable  aeeuraey  by  only  evaluating 
J{A)  and  z{A,R)  at  sueh  points.  We  modify  Proposition  2  slightly  so  that  it  eorresponds  to 
evaluating  y(zl)  only  at  values  A  =  Uk. 

Proposition  3:  Given  a  fixed  A  =  Uk^ox  some  k,  we  ean  solve  MXFI(i?),  for  some 
as  yet  unspeeified  R,  through  LRD(yl,7?),  as  follows: 

1.  Find  a  minimum-eapaeity  Ac  cut  in  G  with  arc  capacities  xxim{uj^,Ar^}\/k  G  A 
where  rk  =  \  for  all  interdietable  ares  k,  and  rk  =  otherwise. 
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2.  Define  the  interdietion  set,  Aj  c  Ac,  where  k  &  A  ^  if  and  only  if  A:  is  a  member 

of  Ac  sueh  that  . 

3.  Let  R  =  \  Ai\  and  define  Xk{R)  =  1  for  all  A:  g  ^4/  and  Xk{R)  =  0,  otherwise.  x{R) 
solves  MXFI(i?)  and  the  slope  of  z{A,R)  is  zero  for  all  A&  {u^.  -£,Ui^)  for  some  £  >  0. 

(We  might  aetually  need  to  perturb  are  eapaeities  to  ensure  that  sueh  an  £  exists.) 

4.  If  yl  =  for  some  k'  G  Ai,  then  let  Ai  =  Aj-  {k'},  R  =  R-  1,  and  define  Xk(R)  =  1 

for  all  A:  G  Aj  and  x^R)  =  0,  otherwise.  x{R)  solves  MXFI(i?)  and  the  slope  of  z{A,R)  is 
zero  for  all  A&  +  e)  for  some  £  >  0  (although  are  eapaeities  might  need  to  be 

perturbed  to  ensure  there  exists  sueh  an  £).  | 


This  suggest  the  following  outline  of  an  algorithm  for  evaluating  x{K)  for 
different  values  of  R: 

1 .  Order  the  Uk,  u\<  U2<.. .  <  u\a\.  These  are  the  potential  values  for  A  (albeit  not  all- 
inelusive). 

2.  Initialize  R  =  °°. 

3.  For  A:=  1,  ...  ,\A\ 

a.  Let  A=  Uk  and  solve  the  eorresponding  max-flow  problem  to  evaluate y(>^), 
Ai,  R,  and  x{R)  as  in  steps  1-3  of  Proposition  3. 

b.  If  R'  R,  print  R  and  x{R)  and  let  R'  =  R  . 

e.  If  yl  =  for  some  k'  e  Ai,  then 

i.  Let  Aj  =  Aj  -  {k'},  R  =  R-  I,  and  define  x{R)  aeeordingly,  all  as  in 
step  4  of  Proposition  3. 

ii.  Print  R  and  x{R)  and  let  R'  =  R  . 
d.  If  i?  =  0,  halt. 
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The  algorithm  above  first  finds  a  solution  x(ri)  to  MXFI(ri)  where  r\  is  the 
number  of  ares  in  a  minimum  cardinality  cut  consisting  of  only  interdictable  arcs:  If  the 
interdictor  had  r\  units  of  interdiction  resource,  he  could  reduce  the  post-interdiction 
maximum  flow  to  zero.  Then,  for  the  same  X  we  may  find  a  solution  to  MXFI(ri-l)  as  in 
step  4  of  Proposition  3.  If  not,  as  X  increases  to  ui  or  some  larger  value,  a  solution  x(r2) 
to  MXFI(r2)  is  found  for  some  r2<  r\.  Unfortunately,  it  may  be  that  r2  <  ri  -  1  and  the 
algorithm  will  have  failed  to  solve  MXFI(ri-l).  Similarly,  the  algorithm  may  solve 
MXFI(r/!)  but  fail  to  solve  MXFI(r/j-i)  and  maybe  MXFI(r/j-2)  and  so  on.  This  issue  will 
be  investigated  further,  but  first  we  show  how  multiple  arcs  with  the  same  capacity  are 
handled. 

C.  SCALING  AND  PERTURBING  ARC  CAPACITIES 

We  have  earlier  assumed  that  all  Uk  are  unique.  If  they  are  not  unique,  we  may 
encounter  a  network  without  easily  identifiable  interdiction  sets  as  Derbes  discovered  in 
his  work.  Derbes’  basic  method  fails  to  find  an  optimal  solution  to  MXFI(i?)  when  it 
cannot  identify  a  set  of  arcs  to  interdict  that  consumes  the  entire  interdiction  resource 
budget  exactly.  One  such  failure  can  occur  when  optimal  cuts  found  in  LRD(/l,i?)  are 
composed  of  a  number  of  arcs  with  equal  capacity.  For  example,  the  network  in  Figure  3 
consists  of  seven  interdictable  arcs  in  parallel  between  the  source  and  the  sink  with  rk  =  \ 
and  Uk=  u  Vk.  Suppose  R  =  4.  There  is  only  one  cut,  but  for  X  >  0,  the  solution  to 
MXFI(i?)  from  Proposition  1  interdicts  all  arcs,  which  is  infeasible,  or  no  arcs,  rather  than 
the  i?  =  4  arcs  that  are  optimal.  Thus,  Lagrangian  relaxation  will  never  find  an  optimal 
solution  to  this  problem. 

In  the  context  of  our  algorithm  which  attempts  to  solve  MXFI(i?)  for  all  integer 
values  of  i?  in  a  certain  range,  this  problem  leads  to  unnecessary  gaps  in  the  values  of  R 
for  which  a  solution  can  be  obtained.  In  this  example,  our  algorithm  would  solve 
MXFI(i?)  for  i?  =  0  and  R  =  l,  but  not  for  any  value  in  between. 
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Figure  1 .  A  network  without  an  easily  identifiable  interdiction  set  when  =  1  for  all  arcs 

k.  All  labeled  arcs  have  the  same  capacity  u  and  are  interdictable;  all  unlabeled 
arcs  have  infinite  capacity  and  are  uninterdictable. 

To  cope  with  this  problem,  Derbes  (1997)  seales  up  all  eapaeities  by  a  faetor  of 
10^  and  then  adds  a  uniform  random  integer  from  the  range  [1,100].  In  practiee  this 
works  well,  but  it  is  still  possible  to  have  two  ares  with  identieal  perturbed  eapaeities, 
whieh  we  would  like  to  avoid. 

We  deal  with  the  problem  of  multiple  arcs  with  identical  capacities  by  simply 
scaling  up  arc  capacities  and  adding  small,  unique  integer  amounts  to  those  ares  with 
identieal  eapaeities.  This  enables  the  algorithm  is  able  to  differentiate  between  ares  that 
would  otherwise  appear  identieal.  The  algorithm  we  use  to  aeeomplish  this  is 

1 .  Sort  the  arcs  so  that  u\<  M2  <  •  •  •  <  u\a\  ; 

2.  10^  M/tforall  kGyf  ; 

3.  For  k  =  2io\A  \  {If  {uk<  Uk-i  )  let  ut  <—  Uk-i  +  1;  } 

We  store  integer  are  eapaeities  as  “long  int”  integers  in  C,  whieh  have  32  bits  of 
aeeuraey  and  ean  represent  integers  somewhat  larger  than  2.1x10^.  Therefore,  our  sealed 
and  perturbed  eapaeities  will  be  representable  as  long  as  lO^Uk  +  \  A  \  <  2.1x10^.  Making 
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the  reasonable  assumption  that  |  ^4  |  <  10  ,  we  ean  be  assured  that  eapaeities  as  large  as  Uk 
=  (2.1x10^  -10*)/10^  =  2x10"^  are  not  too  large:  Our  largest  eapaeities  are  49.  If 
eapaeities  larger  than  2x10"^  are  desired,  the  long  int  integers  ean  be  replaced  by  64-bit 
“long  long  int”  integers.  (64-bit  arithmetic  will  be  carried  out  correctly,  directly  in 
hardware  or  indirectly  through  software,  depending  on  the  computer  used.) 

This  scaled-perturbation  technique  will  work  fine  as  long  the  original  arc 
capacities  are  not  too  large — this  has  just  been  established  for  our  data — and  as  long  as 
there  are  not  “too  many”  arcs  with  identical  capacities.  In  particular,  the  technique  will 
not  cause  an  error  here  unless  for  two  cuts  and  A^^ ,  S  <  E  Mj. ,  but 

kG  kG  ^2 

i.e.,  the  original  capacity  of  cut  A^,^  is  less  than  the  capacity  of  cut  A^,^ , 

kGA^^  kGA^^ 

but  using  the  scaled  and  perturbed  capacities  u[  reverses  the  relationship.  This  cannot 
happen  if  the  total  perturbation  added  to  any  minimal  cut  Ac  does  not  exceed  the  scale 
factor  of  10^.  An  upper  bound  on  the  total  perturbation  added  to  Ac  is  m\  Ac  \  where  m  is 
the  maximum  number  of  arcs  with  identical  capacities.  Conservative  assumptions  for  our 
data  are  that|  I  <  60  and  m  <  600,  so  this  requirement  is  satisfied. 

After  scaling  and  perturbing  arc  capacities,  the  Lagrangian  algorithm  can  be 
modified  to  identify,  for  a  given  R,  certain  multiple  optimal  solutions.  In  particular,  the 
algorithm  might  encounter  and  identify  interdiction  sets,  which  have  identical  original 
capacities  but  different  perturbed  capacities.  This  would  be  accomplished  in  the  outlined 
algorithm  by  replacing  the  “if  statemenf’  in  3.b  with  a  statement  that  compares  the 
elements  of  the  previous  set  Ai  with  the  current  one,  and  prints  out  the  solution  if  the  two 
sets  are  different.  For  the  sake  of  simplicity,  we  ignore  this  possibility. 


D,  A  DETAILED  ALGORITHM 

In  this  section,  we  present  a  detailed  version  of  the  algorithm  outlined  in  section 
3.B.  This  algorithm.  Algorithm  1,  attempts  to  solve  MXFI(i?)  for  all  integers 
Rg  [0,  I  Amin  |],  where  Amin  is  a  minimum-cardinality  cut  consisting  of  only  interdictable 
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arcs.  But  as  stated,  it  only  solves  for  “normal”  values  of  R,  i.e.,  those  values  of  R  that  are 
not  problematie. 

Algorithm  1  requires  the  solution  of  a  sequenee  of  maximum-flow  problems,  and 
beeause  eapaeities  are  non-deereasing  in  this  sequenee,  the  maximum  flow  y'  in  iteration 
h  is  feasible  for  iteration  h+l.  Thus,  it  makes  sense  effieieney-wise  to  use  a  flow- 
augmenting-path  max-flow  algorithm  that  maximizes  flow  starting  with  the  previous 
iteration’s  maximum  flow;  and  our  implementation  of  Algorithm  1  does  ineorporate  this 
feature.  However,  for  the  sake  of  simplieity,  we  do  not  show  this  in  our  pseudo-eode. 
This  maximum-flow  algorithm  is  represented  by: 


Procedure  FindMaxFlow  (G,  s,  t,  u) 

Input:  Network  G  =  (N^),  souree  node  5,  sink  node  t 
u  integer  are  eapaeities  >  0  \lk^  A 

Output:  a  minimum-eapaeity  eut  cz  A 

y  veetor  of  maximum  are  flows,  >  OVk  G  A 

f  maximum-flow  value 

{ 

This  proeedure  finds  a  max  flow  y,  max-flow  value/  and  identifies  a  minimum- 
eapaeity  eut  A^  using  a  variant  of  the  Edmonds  and  Karp  (1972)  max-flow 
algorithm  as  implemented  by  Derbes  (1997). 

return  ( ,  y,/) 

} 

Algorithm  1  ean  now  be  stated.  Note  that  solutions  are  only  given  in  terms  of  the 
value  for  R,  the  optimal  interdietion  set  Aj  given  R,  and  the  optimal  objeetive  value  to 
MXFI(7?)  whieh  is  z(R). 

Algorithm  1:  Lagrangian-Relaxation  Algorithm  for  MXFI 

Input:  Network  G  =  (N^A),  souree  node  5,  sink  node  t 

u  integer  are  eapaeities  A,  dX\  u^.  assumed  unique 

=  1  for  all  interdietable  ares  Vk  G  H =  0°  otherwise 

Output:  Ai  optimal  interdietion  set  for  eaeh  normal  value  of  7?  G  [0,|  Amin  |]; 

z{R)  optimum  objeetive  value  to  MXFI(7?)  for  eaeh  normal  value  of  R; 

Step  0:  /*  define  and  order  the  eandidate  set  of  A  values  */ 
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^ \fk  &  A', 

Order  the  values  At  such  that  Ai  <A2<  •  •  •  <  A\a\; 

Step  1:  Initialize:  R<^\A\  +l,i?'<— |  A  \  +1,/  <—  0,/'<—  0,yl'<—  0,  /  <—  1 ,  y'<—  0; 

Step  2: 

While  {1<\A  \  and  R^0){ 

A  < —  Ai', 

u[  <—  min{M^,ylr^}  for  all  A:  G  ^4 ;  /*  adjust  arc  capacities  */ 

{A^,y,f)  <—  FindMaxFlow  (G,  s,  t,  u'); 

A/  <—  {k  G  Ac  I  =  A }; 

R<r-  \Al\', 

z  <— /-  AR', 

/*  Uncomment  the  next  block  to  handle  problematic  i? -values  */ 

/* 

if  (i?  -i?'>  1 ){ 

Call  SolveForProblematic(i?,  A,f,  R',  A',f',  A'^); 

} 

*  / 

print  (“  The  optimal  solution  to  MXFI(i?)  for  R  =  ”,  R,  “follows:  ”); 
print  (The  optimal  interdiction  set  and  objective  value  are”,  Aj,  z)  ; 
if  (  Uj^'=  A  for  some  k'  g  Aj){ 

Ai<^  Aj—  {k'}', 

R^R-\; 
z  <—  z  +  ; 

print  (“  The  optimal  solution  to  MXFI(i?)  for  R  =  ”,  R,  “follows:  ”); 
print  (The  optimal  interdiction  set  and  objective  value  are”,  Aj,  z)  ; 

} 

R  R,  A  A  ,  y  <—  y^,  Aj  <—  Aj  ,/'<—/  z'  <—  z; 

/^/  +  1; 


If  the  change  in  i?-values  from  one  iteration  to  the  next,  is  never  larger  than  one. 
Algorithm  1  must  solve  MXFI(i?)  for  all  values  of  R  in  the  range  [0,  |  Amin  |  ]  as  described 
in  section  3.B.  When  the  change  exceeds  one,  there  will  be  one  or  more  missing 
solutions  and  we  will  modify  the  algorithm  by  adding  the  procedure 
SolveForProblematic(),  to  deal  with  this.  This  heuristic  is  described  in  the  next  section. 


E.  PROBLEMATIC  VALUES  OF  R 

We  don’t  know  when  problematic  i?-values  will  arise,  but  we  know  empirically 
that  they  are  likely  to  occur  and  must  therefore  devise  a  method  to  cope  with  them.  We 
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devise  a  teehnique  here  to  evaluate  bounds  and  feasible  interdietion  sets  for  the 
problematie  values  of  R  by  ealling  the  proeedure  SolveForProblematie()  whieh  is 
eommented  out  in  Algorithm  1 .  We  eall  the  full  version  of  the  algorithm  Algorithm  2. 

When  Algorithm  2  finds  a  solution  for  i?  =  r  and  R  =  r  -  2,  but  not  for  i?  =  r  -  1 ,  it 
ealls  SolveForProblematieO  to  devise  a  solution  MXFI(r-l):  It  “uninterdiets”  the  are 
with  the  least  eapaeity  in  the  “i?=r  solution”.  This  yields  an  upper  bound  UB  >  z(r-l),  of 
eourse,  and  we  eompute  a  lower  bound  T5  =  max{j{Ar-2)  -  Ar-2{r-\),  J{Ar)  -  Ar{r-l)  } 
where  A.r  -  2  maximizes  z{A,r-2)  and  Ar  maximizes  z{A,r)  (and  led,  respeetively,  to  the 
solutions  of  MXFI(r-2)  and  MXFI(r)).  Absolute  and  relative  optimality  gaps  are 
eomputed  for  the  solution  as  AbsGap  =  UB  -  LB  and  RelGap  =  {lOOViAixAbsGapILB. 

In  praetiee,  Algorithm  2  also  ereates  a  heuristie  solution  fori?  =  r-1  by  adding  an 
interdietion  to  the  “i?=r-2  solution.”  We  do  not  show  this  in  the  statement  of  the 
algorithm  for  the  sake  of  simplieity.  For  the  sake  of  eomputational  effieieney,  we  have 
not  attempted  to  optimize  the  lower  bound  LB,  i.e.,  we  have  not  tried  to  maximize 
z{A,r-\).  To  do  this  would  require  that  we  explore  the  region  {Ar  -  2,  Ar)  more  fully  and 
solve  more  maximum- flow  problems. 

When  a  sequenee  of  two  or  more  problematie  i?-values  is  eneountered,  we  simply 
apply  the  above  proeedure  repeatedly;  this  will  be  shown  by  the  seeond  example  below. 
In  these  examples,  the  notation  AAr)  denotes  the  interdietion  set  identified  for  R  =  r. 
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Figure  2.  Sample  Network,  NET5,  to  Demonstrate  Problematic  7?-values.  Numbers  on 
arcs  are  capacities  and  letters  are  labels.  Unlabeled  arcs  cannot  be  interdicted, 
i.e.,  rt  =  All  labeled  arcs  have  r*  =  1. 

To  illustrate  SolveForProblematicQ,  consider  NET5  shown  in  Figure  2. 
Algorithm  2  finds  a  solution  for  i?  =  9  and  then  for  R  =  1,  but  not  for  R  =  S.  For 
reference:  z(9)  =  100;  =  2703,  where  Aj  =  300  optimizes  z{A,7)  =  z(7);  the  R=9 

solution  is  A](9}  =  {a,b,c,d,e,f,g,h,k};  and^zlg)  =1918  where  Ag  =  202.  To  solve  MXFI(8) 
approximately,  we  carry  out  the  following  steps: 

i.  Ai^)=Ai9)  -  /cmin=  {a,b,c,d,e,f,g,h,k}-  {k}  =  {a,b,c,d,e,f,g,h}, 

ii.  UB  =  z(9)  +  Mm  =  100  +  300  =  400, 

iii.  LB  =  max{/i7)  -  Ai^,f{Ag)  -  AgS}  =  max {2703-300(8),  1918-202(8)}  =  303, 
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iv.  AbsGap  =  UB  -  LB  =  400  -  303  =  97,  and 

V.  RelGap  =  \0()yAbsGAPILB  =  32. 

The  problem  of  Figure  2  also  provides  an  example  of  two  eonseeutive 
problematie  7?-values.  In  partieular,  Algorithm  2  solves  exaetly  for  7?  =  7  (/I?  =  502)  and 
for  7?  =  4  (/I4  =  600),  but  not  for  7?  =  6  and  R  =  5.  For  referenee,  z(7)  =  603;  =  4117, 

f{A4)  =  4604,  where  yl4  =  600  optimizes  z(yl,4)  =  z(4);  and  A]{7)  =  {m,n,o,q,r,u,x} .  Then, 
the  algorithm  will  eompute: 

1.  For  7?  =  6, 

i.  Aj(6)  =  {m,n,o,q,r,u,x}  -  =  {m,n,o,q,r,u,x}  -  {m}  =  {n,o,q,r,u,x}, 

ii.  UB  =  z(7)  +  Mm  =  603  +  601  =  1204, 

hi.  LB  =  max{/(i4)  -  ^a6,A^i)  -  >^76}  =  {4604  -  600(6),  4117  -  502(6)}  =  1105, 

iv.  AbsGap  =  UB  -  LB  =  1204  -  1 105  =  99,  and 

V.  RelGap  =  lOOxAbsGAPILB  ='^.9. 

2.  For  7?  =  5, 

i.  Ai=  {n,o,q,r,u,x}  -  ^min  =  {n,o,q,r,u,x}  -  {u}  =  {n,o,q,r,x}, 

ii.  6/5  =  z(6)  + Mu  =  1204  +  602  =  1806, 

hi.  LB  =  max{/(i4)  -  -  ^1^)  {4604  -  600(5),  4117  -  502(5)}  =  1607, 

iv.  AbsGap  =  UB  -  LB  =  1806  -  1607=  199,  and 

V.  RelGap  =  lOOxAbsGAPILB  =  12.4. 

This  example  also  illustrates  that  our  methodology  does  not  neeessarhy  eompute 
the  best  lower  bound.  For  instanee,  the  lower  bound  we  obtain  for  7?  =  5  is  1607,  but  the 
best  lower  bound  is  max^  z(yl,5)  =  1670  whieh  oeeurs  at  h.  =  535.  Algorithm  2  will  never 
diseover  this  beeause  it  only  evaluates  z{A,R)  for  A  equal  to  some  are  eapaeity;  and  535  is 
not  an  are  eapaeity  in  this  example. 

SolveForProblematie( )  is  now  speeified  in  detail  here: 
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Procedure  SolveForProblematic  {R,  A,f,  R',  A',f',  A',) 

Input:  R  current  i?-value  correctly  solved 

A  value  of  parameter  that  maximized  z{A,R)  to  solve  MXFI(i?) 

/  maximum  flow^zl) 

R'  previous  i?-value  correetly  solved  for,  {R'  >  R  +  1) 

X  value  of  parameter  that  maximized  z{A,R')  to  solve  MXFI(i?') 

Output:  A'j  interdiction  sets  that  are  approximate  solutions  for  problematic 

i?- values  in  the  range  R'-l,  R'-2,  ,  i?+l 

UB  upper  bound  for  problematic  value  of  R 

LB  lower  bound  value  for  problematic  value  R 

AbsGAP  absolute  optimality  gap  UB  -  LB 
RelGAP  relative  optimality  gap  ((100%)x(175  -  LB)ILB 

Step  0: 

UB  ^  z; 

Step  1: 

For  (  r  =  i?'  -  1  down  to  R  +  1  )  { 
argminw^ ; 

keAj 

LB  <— max{/-  Ar,f'-  Xr  }; 

UB  < —  UB  +  j 

"■min 

AbsGAP  ^UB-  LB; 

RelGAP  ^  im-xAbsGAPILB; 

^'l  ^  ^min  ’ 

print  (“  For  problematic  value  R  =  ”,  r,  “the  approximate  soln.  is:  ”); 
print  ( 4  ,UB,  LB,  AbsGAP,  RelGAP  ); 

) 

E.  COMPLEXITY  OF  THE  ALGORITHM 

Finding  the  minimum-capaeity  cut  after  finding  the  maximum-flow  requires 
0(1^4 1)  time  at  worst.  The  sort  of  the  values  for  Ak  has  complexity  O(mlogm).  If  no 
problematie  values  of  R  are  observed,  the  work  in  the  algorithm  is  elearly  dominated  by 
the  need  to  solve  perhaps  as  many  as  m  maximum- flow  problems.  In  this  case,  the 
complexity  of  the  overall  algorithm  is  0{mg{m,n))  where  the  max-flow  algorithm  has 
complexity  0{g{m,n)).  At  most  m  problematic  values  of  R  can  be  encountered  and  the 
work  required  to  deal  with  these  is  at  most  0{m)  each.  Thus,  the  overall  eomplexity  of 
dealing  with  problematic  i?-values  is  0(m  ).  The  overall  complexity  remains  0(mg(m,n)) 
since  0{g{m,n))  is  eertainly  worse  than  0{m).  In  practice,  we  may  achieve  better 
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performance  than  this  because  we  are  solving  a  sequence  of  related  maximum  flow 
problems,  and  not  solving  each  problem  from  scratch. 
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IV.  COMPUTATIONAL  RESULTS 


We  describe  the  networks  tested  in  the  first  part  of  this  chapter  and  present 
computational  results  in  the  second  part. 


A,  TEST  NETWORK  DESIGN 

We  have  tested  our  algorithm  using  five  artificial,  directed  networks;  see  Table  1. 
The  structure  of  these  networks  should  be  similar  to  the  structure  of  road  networks  of 
interest.  These  networks  are  based  on  x  m  grid  networks  where  n\  is  the  number  of 
nodes  on  the  horizontal  axis  and  n2  is  the  number  of  nodes  on  the  vertical  axis. 
Horizontal  arcs  are  oriented  from  “west  to  east’  and  all  nodes  directly  above  or  below 
another  are  connected  with  a  pair  of  anti-parallel  arcs.  Anti-parallel  arcs  are  separately 
interdictable.  We  also  have  arcs  oriented  in  the  southeast  and  northeast  directions.  All 
westernmost  nodes  are  connected  to  a  source  node  with  artificial,  uninterdictable, 
infinite-capacity  arcs  and  all  easternmost  nodes  are  connected  to  a  sink  node  with 
artificial,  uninterdictable,  infinite-capacity  arcs.  The  arc  capacities  Uk  are  randomly 
drawn  from  the  discrete  uniform  distribution  on  [1,49]  and  r^=\  is  assumed  for  all 
interdictable  arcs  k. 

We  base  our  implementation  of  Algorithm  2  on  Derbes  (1997),  which  uses  an 
improved  version  of  the  Edmonds  and  Karp  (1972)  maximum- flow  algorithm.  His  code 
only  attempts  to  solve  MXFI(7?)  for  single  values  for  7?;  ours  attempts  to  solve  for  all 
values  7?  G  [0,  |  Hmin  |]  where  Hmin  is  a  minimum-cardinality  cut  in  the  network  consisting 
of  only  interdictable  arcs. 
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Network 

«i 

«2 

n 

m 

/* 

NET27 

5 

5 

27 

86 

238 

NET51 

7 

7 

51 

188 

361 

NETlOl 

11 

9 

101 

412 

414 

NET  122 

8 

15 

122 

499 

765 

NET402 

20 

20 

402 

1826 

1080 

Table  1.  Statistics  for  Test  Networks,  n  is  the  total  number  of  nodes  and  m  the  total 
number  of  nodes.  The  networks  are  based  on  ni  x  112  grids  of  nodes.  The 
number  of  interdictable  arcs  is  m  -  2n.  J*  is  the  uninterdicted  maximum  flow, 
where  arc  capacities  are  uniformly  distributed  random  integers  in  the  range 
[1,49]. 

B,  RESULTS 

Test  results  are  summarized  in  Table  2,  whieh  explieitly  shows  only  the  results  for 
problematie  i?-values  sinee  MXFI(i?)  is  solved  exaetly  for  other,  normal  values.  We  are 
primarily  interested  the  pereentage  of  problematie  i? -values  and  the  aeeuraey  of  results 
for  sueh  values.  The  maximum  flow  value  without  interdietion  is  indieated  hy  f*  in  the 
tables.  The  maximized  value  of  the  Lagrangian  funetion  is  z{A,R)  while  z{A,R)  is  the 
upper  bound  value,  whieh  equals  z{A,R)  when  R  is  not  problematie  but  is  obtained  from 
the  heuristie  when  it  is. 

Table  2  shows  that  no  problematie  i?-values  oeeur  for  the  first  two  problems  and 
that,  in  general,  problematie  values  oeeur  less  than  25%  of  the  time.  Furthermore, 
absolute  optimality  gaps  are  always  small,  being  at  most  3%  of  the  uninterdieted 
maximum  flow.  NET402  exhibits  10  out  of  58  i? -values  as  problematie  all  of  whieh  have 
relative  optimality  gaps  of  less  than  2.7%.  Large  relative  optimality  gaps,  33.3%  and 
20.0%,  oeeur  in  two  instanees  for  NET122,  but  these  also  depend  on  the  small  value  of 
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the  lower  bound.  The  rest  of  the  problematie  i?-values  for  this  network  have  relative  gaps 
of  less  than  1.1%,  whieh  seem  very  reasonable. 

In  Tables  3  and  4,  we  explore  ehanges  in  optimality  gaps  for  one  of  the  test 
networks,  NET  122,  as  a  funetion  of  the  eoarseness  of  the  are-eapaeity  distribution.  When 
the  eapaeity  distribution  is  eoarse,  i.e.,  there  are  few  are  eapaeities,  the  eurves  for  z{A,R) 
will  not  be  evaluated  at  many  points  and  we  ean  expeet  that  lower  bounds  will  not  be 
very  aeeurate.  Furthermore,  fewer  potential  solutions  will  be  seen  and  thus  we  expeet 
upper  bounds  to  be  worse.  We  note  that  the  largest  problem,  NET402,  requires  384 
seeonds  to  solve  if  eaeh  maximum-flow  problem  is  solved  from  serateh,  so  solving  the 
maximum- flow  problem  in  iteration  h+\  starting  with  the  flow  from  iteration  h,  as  we  do, 
ean  yield  a  substantial  improvement  in  run  time. 

In  Table  3,  NET122  is  evaluated  with  eapaeities  randomly  drawn  from  {16,32}, 
{8,16,24,32},  and  {4,8,12,16,20,24,28,32}  (denoted  16[1,2],  8[1,4],  4[1,8],  respeetively); 
in  Table  4  it  is  evaluated  with  eapaeities  drawn  from  {2,4, 6..., 32},  and  {1,2,..., 32} 
(denoted  2[1,16],  [1,32],  respeetively).  As  suspeeted,  these  tables  indieate  that  more 
eoarsely  distributed  are  eapaeities  lead  to  poorer  results:  The  eoarsest  distribution  has  the 
largest  absolute  and  relative  errors  and  the  least  eoarse  the  smallest.  It  also  appears  from 
the  results  that  when  there  is  a  non-zero  optimality  gap,  a  lower  bound  on  this  gap  is  the 
smallest  non-zero  differenee  between  are  eapaeities.  In  partieular,  the  smallest  non-zero 
gap  for  16[1,2]  is  16,  for  8[1,4]  is  8,  and  so  on.  Also,  we  see  that  gaps  appear  to  be 
monotonieally  non-deereasing  in  a  eontiguous  sequenee  of  problematie  i?-values.  For 
instanee,  eonsider  the  sequenee  of  R-values  from  19  down  to  10  in  16[1,2]  and  similar 
sequenees  in  the  other  problems. 
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Network 

R 

Run  Time 
fcDU  sec.l 

Problematic 

R  -Values 

z(lR) 

z(A,R) 

Abs. 

Gap 

Rel. 

Gap 

NET27 

238 

[1,131 

0.1 

none 

NA 

NA 

NA 

NA 

NET51 

361 

[1,191 

0.5 

none 

NA 

NA 

NA 

NA 

NETlOl 

414 

[1,251 

1.0 

13 

75 

73 

2 

2.7 

NET  122 

765 

[1,431 

1.0 

41 

4 

3 

1 

33.3 

40 

6 

5 

1 

20.0 

23 

166 

165 

1 

0.6 

22 

184 

182 

2 

1.1 

21 

200 

199 

1 

0.5 

17 

277 

276 

1 

0.4 

16 

299 

297 

2 

0.7 

15 

320 

319 

1 

0.3 

11 

410 

408 

2 

0.5 

10 

432 

431 

1 

0.2 

NET402 

1080 

[1,581 

16.0 

47 

37 

36 

1 

2.7 

46 

42 

41 

1 

2.4 

38 

113 

112 

1 

0.9 

33 

175 

173 

2 

1.1 

32 

190 

186 

4 

2.1 

31 

202 

200 

3 

1.0 

22 

350 

349 

1 

0.3 

21 

373 

371 

2 

0.5 

20 

396 

393 

3 

0.8 

19 

417 

416 

1 

0.2 

Table  2.  Summarized  Test  Results  for  the  Five  Test  Networks.  /*  is  the  uninterdieted 
value  for  maximum  flow.  “NA”  indicates  “not  applicable,”  used  because 
MXFI(R)  solves  exactly  for  all  values  off?  for  NET27. 
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R 

Run  Time 
tcDU  sec.l 

Problematic 

R  -Values 

z(A.R) 

z(A,R) 

Abs. 

Gap 

Rel. 

Gap 

16[1,2]  912 

[1,43] 

1.5 

19 

400 

384 

16 

4.2 

18 

432 

400 

32 

8.0 

17 

464 

416 

48 

11.5 

16 

496 

432 

64 

14.8 

15 

528 

448 

80 

17.8 

14 

544 

464 

80 

17.2 

13 

560 

496 

64 

12.9 

12 

576 

528 

48 

9.9 

11 

592 

560 

32 

5.7 

10 

608 

592 

16 

2.7 

8  [1,4]  672 

[1,43] 

1.5 

28 

128 

120 

8 

6.7 

27 

136 

128 

8 

6.3 

18 

280 

272 

8 

2.9 

17 

304 

288 

16 

5.5 

16 

328 

304 

24 

7.9 

15 

344 

320 

24 

7.5 

14 

360 

336 

24 

7.1 

13 

376 

352 

24 

6.8 

12 

392 

368 

24 

6.5 

11 

408 

384 

24 

6.3 

10 

424 

408 

16 

3.9 

9 

440 

432 

8 

1.8 

4[1,8]  584 

[1,43] 

1.5 

21 

172 

168 

4 

2.4 

14 

284 

280 

4 

1.4 

13 

300 

296 

4 

1.4 

12 

316 

312 

4 

1.3 

11 

332 

328 

4 

1.2 

10 

348 

344 

4 

1.2 

9 

364 

360 

4 

1.1 

3 

496 

492 

4 

0.8 

Table  3.  Partial  Test  Results  for  NET122  with  Restrieted  Arc  Capacities.  16[1,2] 

indicates  capacities  uniformly  distributed  from  the  set  {16,32},  8[1,4]  from  the 
set  {8,16,24,32}  and4[l,8]  from  the  set  {4,8,.  ..,32}. 


Wk 

R 

Run  Time 
fcDU  sec.) 

Problematic 

R-Values 

z(A.K) 

z{A,R) 

Abs. 

Gap 

Rel. 

Gap 

2[1,16] 

530 

[1,43] 

1.5 

21 

144 

142 

2 

1.4 

20 

156 

154 

2 

1.3 

16 

212 

210 

2 

1.0 

15 

226 

224 

2 

0.9 

14 

240 

238 

2 

0.8 

11 

288 

286 

2 

0.7 

[1,32] 

503 

[1,43] 

1.5 

22 

121 

120 

1 

0.8 

21 

133 

132 

1 

0.8 

16 

199 

197 

1 

1.0 

15 

213 

211 

2 

0.9 

14 

226 

225 

1 

0.4 

13 

240 

239 

1 

0.4 

12 

254 

253 

1 

0.4 

11 

268 

267 

1 

0.4 

Table  4.  Partial  Test  Results  for  NET122  with  Restrieted  Arc  Capacities.  2[1,16] 

indicates  capacities  uniformly  distributed  from  the  set  {2, 4, ...,32}  and  [1,32] 
from  the  set  {1,2,.  ..,32}. 
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V.  CONCLUSIONS  AND  FUTURE  WORK 


A,  SUMMARY 

This  thesis  has  studied  a  speeial  solution  teehnique  for  solving  or  approximately 
solving  the  maximum-flow  network-interdietion  problem  (MXFI).  In  MXFI,  a  network 
user  attempts  to  maximize  flow  of  a  eommodity  through  the  network,  while  an  interdietor 
destroys  network  ares,  using  limited  interdietion  resourees,  to  minimize  this  maximum 
flow. 

MXFI  ean  be  modeled  as  a  binary  integer  problem  and  solved  but  this  approaeh 
ean  be  slow.  A  Lagrangian-relaxation  teehnique  developed  by  Derbes  (1997)  has  shown 
promise  for  solving  the  problem  more  quiekly,  and  this  thesis  has  extended  Derbes’ 
teehnique  for  solving  MXFI  for  all  values  of  i?  in  a  speeified  range.  It  is  assumed  that 
one  unit  of  resouree  is  required  to  interdiet  any  interdietable  are,  so  we  are  eoneemed 
only  with  integer  values  of  R. 

In  effeet,  our  basie  algorithm,  Algorithm  1,  determines  the  breakpoints  in  the 
pieeewise-linear,  eoneave  Lagrangian  funetion,  whieh  ean  only  oeeur  when  the 
Lagrangian  multiplier  A  equals  the  eapaeity  of  some  are.  For  eaeh  linear  pieee  of  the 
funetion,  MXFI  is  solved  exaetly  for  some  value  of  R,  with  inereasing  values  of  A 
eorresponding  to  smaller  values  of  R.  But  the  algorithm  may  leave  gaps  in  the  sequenee 
of  7?-values  for  whieh  the  problem  is  solved.  Algorithm  2  eombines  Algorithm  1  and  a 
heuristie  to  fdl  in  these  gaps,  usually  quite  sueeessfully. 

We  have  performed  basie  testing  of  Algorithm  2,  eoded  in  C,  using  five  test 
networks  ranging  in  size  from  27  nodes  and  86  ares  to  402  nodes  and  1826  ares.  Are 
eapaeities  are  uniformly  distributed  random  integers  in  the  range  [1,49].  All  tests  are 
performed  on  a  700  MHz  Pentium  III  personal  eomputer  with  The  Mierosoft  Windows 
Millennium  operating  system  and  Mierosoft  Visual  C++  eompiler.  The  largest  network 
is  solved  in  16  seeonds.  Absolute  optimality  gaps  are  always  small,  at  most  two  pereent 
of  the  uninterdieted  maximum  flow.  Relative  gaps  are  typieally  small,  a  few  pereent,  but 
ean  be  large  if  the  post-interdietion  maximum  flow  is  small. 
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Additional  testing  evaluated  one  of  the  networks  with  eapaeities  distributed  more 
and  less  eoarsely.  In  partieular,  eapaeities  are  drawn  randomly  from  {16,32}, 
{8,16,24,32},. ..,32}.  We  eonelude  that  absolute  and  relative  optimality  gaps  tend 
to  be  larger  with  the  eoarser  eapaeity  distributions. 

B,  FUTURE  WORK 

Substantial  effieieney  is  gained  in  our  solution  proeedure  by  solving  the 
maximum- flow  problem  in  iteration  h+\  starting  with  the  solution  from  iteration  h;  flow- 
augmenting-path  maximum-flow  algorithms  work  well  in  this  eontext.  However,  there 
exist  more  effieient  flow-augmenting-path  algorithms  than  ours,  and  it  might  be  worth 
investigating  sueh  algorithms  for  the  purpose  of  improving  overall  run  times. 

Further  researeh  on  this  problem  may  foeus  on  the  solution  methods  of 
problematie  i?-values.  We  have  devised  a  heuristie  to  deal  with  these,  but  an  exaet 
braneh-and-bound  algorithm  eould  exploit  the  lower  bounds  provided  by  our  Lagrangian 
algorithm.  We  also  know  that  our  lower  bounds  are  not  always  maximal,  and  this  issue 
needs  to  be  investigated  further:  Improved  bounds  eould  be  obtained  by  simply  defining 
a  finer  grid  for  A  than  the  one  we  use  (whieh  eorresponds  to  are  eapaeities),  but  of  course 
this  would  add  to  the  computational  burden.  The  finer  grid  might  also  lead  to  better 
feasible  solutions,  too.  It  may  also  be  possible  to  reformulate  the  problem  so  that 
Lagrangian  relaxation  yields  smaller  optimality  gaps.  Of  course,  the  technique  should  be 
generalized  so  that  the  resource  required  to  interdict  an  arc,  rk,  is  a  general  integer  rather 
than  restricting  =  1 . 

The  model  MXFI  can  be  extended  to  dynamic  networks  by  adding  traversal  time 
for  each  arc  in  the  networks.  The  model  can  also  be  extended  to  consider  stochastic  arc 
capacities  and/or  uncertain  interdiction  successes  (Cormican,  Morton  and  Wood  1995). 
Perhaps  our  Lagrangian  technique  could  be  used  in  the  solution  procedure  for  these 
problems. 

Golden  (1978),  Israeli  (1999),  Israeli  and  Wood  (1999)  and  Wevley  (1999)  study 
another  category  of  network  interdiction  model  where  an  interdictor  attempts  to  interdict 
(destroy  or  lengthen)  arcs,  using  limited  interdiction  assets,  to  maximize  the  length  of  a 
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shortest  s-t  path.  While  these  models  are  similar  to  MXFI,  that  earlier  work  is  not 
relevant  to  this  thesis.  However,  the  Lagrangian-relaxation  teehnique  developed  here 
may  have  an  analog  for  shortest-path  interdietion  and  this  may  be  worth  investigating. 
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